xend = 0.1
xstart=0;

a=zeros(1,6);

T=[ 1 0 0 0 0 0 ;
    0 1 0 0 0 0 ;
    0 0 2 0 0 0 ;
    1 1 1 1 1 1 ;
    1 1 2 3 4 5 ;
    0 1 2 6 12 20 ;    
  ];

a=T^-1*[xstart;0;0;xend;0;0]
